A Radio Pulsar Spinning at 716 Hz 
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We have discovered a 716-Hz eclipsing binary radio pulsar in the globular 
cluster Terzan 5 using the Green Bank Telescope. It is the fastest-spinning 
neutron star ever found, breaking the 23-year-old record held by the 642-Hz 
pulsar B1937+21. The difficulty in detecting this pulsar, due to its very low flux 
density and high eclipse fraction (~40% of the orbit), suggests that even faster- 
spinning neutron stars exist. If the pulsar has a mass less than 2 M , then its 
radius is constrained by the spin rate to be < 16 km. The short period of this 
pulsar also constrains models that suggest gravitational radiation, through an 
r-mode instability, limits the maximum spin frequency of neutron stars. 

While the majority of neutron stars are observed to rotate slower than a few times a second, 
those in binary systems can reach spin rates of hundreds of times a second through the transfer 
of angular momentum from their companion star (f/JEl). Some of these neutron stars, termed 
millisecond pulsars, are persistent radio sources whose emission is modulated at the star's spin 
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frequency. Determining the maximum achievable rotation rate of a neutron star is important 
for a variety of astrophysical problems, ranging from understanding the behaviour of matter 
at supra-nuclear densities, to estimating the importance of neutron stars as gravitational wave 
sources for current and upcoming gravitational wave detectors. For over 23 years, the 642-Hz 
pulsar B 1937+21, the first millisecond pulsar ever found, has been the fastest-spinning neutron 
star known ©. It has been argued that faster ones are exceedingly rare, if they exist at all (4). 

Per unit mass, globular clusters have many more millisecond pulsars than the Galactic disk. 
This is due to the extremely high stellar densities in their cores (10 4 — 10 6 pc~ 3 ), which promote 
the creation of binary systems © where a neutron star is spun-up (or "recycled") to rotate 
hundreds of times a second (2101. We have searched the massive and dense globular cluster 
Terzan 5 for millisecond pulsars using the National Radio Astronomy Observatory's © 100-m 
Green Bank Telescope (GBT). Our searches have thus far uncovered 30 millisecond pulsars in 
Terzan 5 (7), in addition to the three previously known pulsars in this cluster (Iffll9ll70l). The 
discovery of 21 pulsars in Terzan 5 was presented in (f77t. Following that paper, an additional 
nine pulsars have been found in searches of our monitoring observations. These will be reported 
elsewhere. Terzan 5 has the largest known population of millisecond pulsars of any globular 
cluster, roughly one quarter of the entire population of globular cluster millisecond pulsars, and 
the five fastest rotating pulsars in the Galactic globular cluster system (0. Among the newest 
discoveries in Terzan 5 is PSR J 1748— 2446ad, a 716-Hz eclipsing binary millisecond pulsar, 
which is the fastest-spinning neutron star ever found. 

Observations and data analysis. We discovered PSR J1748— 2446ad in 10 November 
2004 observations of Terzan 5 and confirmed it in 8 January 2005 observations using the Pulsar 
Spigot backend (|72")l on the GBT. All observations employed the Spigot with 600 MHz of usable 
bandwidth centered at 1950 MHz, 768 spectral channels, and 81.92-/is sampling. Observations 
were generally 6—7 hours in length and taken at roughly monthly intervals starting June 2004. 
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In addition, there was a more closely spaced set of observations in early May 2005. 

The discovery observation showed that the pulsar is part of a binary system and is eclipsed 
by its companion; both of these properties restricted our ability to detect the pulsar in our 
monitoring sessions. Nonetheless, we have now detected the pulsar in at least 18 out of the 
30 multi-hour observations taken thus far (see Fig 1 for pulse profile). We have derived a 
reliable orbital ephemeris (see Tab 1) by initially modelling the pulse phase delays of a few 
good detections with a simple sine function and then refining the model by fitting pulse times of 
arrival to a simple Keplerian orbital model, with arbitrary pulse phase offsets between observing 
epochs. This ephemeris allowed us to detect the pulsar in many observations where it was not 
initially identified through a periodicity search. 

The pulsar is in a highly circular 26-hr orbit with a >O.14M companion, and is eclipsed 
for ~40% of its orbit at 2 GHz. Such a large eclipse fraction, corresponding to an eclipse re- 
gion with physical size ~5— 6 R Q is extremely rare for such a relatively wide orbit (separation 
between the pulsar and companion of ~4— 5R ). The companion may be a bloated main- 
sequence star, possibly still filling its Roche Lobe, as has been suggested for PSR J1740— 5340, 
a 35-hr binary millisecond pulsar with a >0.21 M Q companion and ~40% eclipse fraction at 
1.4GHz (TOt. The eclipse properties are also similar to those of PSR J1748-2446A, a 1.8- 
hr binary with a 0.089 M minimum mass companion, also located in Terzan 5 ©. Like 
PSR J1740-5340 and PSR J1748— 2446A, there is evidence that the eclipse duration of PSR J1748-2446ad 
is variable, and sometimes lasts longer than 40% of the orbit. On a least two epochs when our 
ephemeris predicts the pulsar should have been visible for several hours, it was not detected 
at all. The pulse signal-to-noise is too low to measure dispersion measure variations on short 
timescales. Future observations should allow a phase-coherent timing solution to be derived, 
which will provide a precise position and observed spin frequency derivative v. Until then, we 
have provided an upper limit on \ v\ (Tab 1). 
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We have verified that this signal is not harmonically related to any of the other known pulsars 
in Terzan 5. We have also carefully investigated the possibility that our searches have identified 
the second harmonic of a new 358-Hz pulsar. When we fold the data at 358 Hz, there are two 
identical peaks (within the resolution and RMS noise level of the Spigot data) separated by 1 80° 
in pulse phase. This is what we expect to see if the pulsar signal is folded at half its intrinsic 
spin frequency, and strongly suggests that 716 Hz is the true spin frequency of the pulsar. The 
results of folding the data at other harmonically related spin frequencies are also consistent 
with the pulsar having a true frequency of 716 Hz. There is also evidence (Fig 1) for a weak, 
but statistically significant interpulse (extra structure in addition to the main pulse peak) when 
the data are folded at 716 Hz. This interpulse, if real, is further evidence that the spin frequency 
is716Hz. 

Lastly, we have simultaneously observed the pulsar using the GASP coherent dedispersion 
pulsar machine (|7?ll75l) . which records only ~l/6 the bandwidth achievable with Spigot but 
which removes all dispersive smearing due to the ionized interstellar medium, resulting in sig- 
nificantly sharper (i.e. narrower) pulse profiles (176b . When the GASP data are folded at 358 Hz, 
two peaks, consistent in shape with each other to within the RMS of the noise, are seen sepa- 
rated by 180° in pulse phase, again indicating that 358 Hz is half the true spin frequency. We 
conclude that PSR J1748-2446ad is indeed a 716-Hz pulsar. 

Implications and discussion. The equation of state of matter at supra-nuclear densities, and 
thus the mass-radius relation for neutron stars, is unknown. If a star is rotating too rapidly for a 
given radius, centrifugal forces will cause it to shed mass. Lattimer & Prakash f/71) derive the 
following equation which, independent of the true equation of state, gives the maximum spin 
frequency for a neutron star with a non-rotating radius R and mass M (assuming this is not close 
to the maximum mass allowed by the equation of state): v max = 1045(M/M Q ) 1 / 2 (10km/R) 3//2 Hz. 
Using this, and assuming a mass less than 2M (which accomodates all measured neutron star 



4 



masses (1771) ) we find an upper limit of 16 km on PSR J1748— 2446ad's radius. We note that this 
constraint applies specifically to PSR J 1748— 2446ad and that slower- spinning pulsars could 
have larger radii. Recently, Li & Steiner (f77fl) have derived a radius range of 11.5— 13.6km 
for a 1 .4 M neutron star, based on terrestrial laboratory measurements of nuclear matter. For 
a 1.4M© neutron star, we find an upper limit of 14.4 km, which is in agreement with their 
result. These radius constraints are more robust than those obtained through observations of 
neutron star thermal emission, which is faint, difficult to measure, and whose characterization 
depends on uncertain atmosphere models (TTPb . Although in principle a radius measurement 
could constrain the unknown equation of state of dense matter, PSR J 1748— 2446ad does not 
rule out any particular existing models, since the pulsar mass is unknown. It is unlikely that a 
mass measurement will be achievable through timing of the pulsar, as the orbit is too circular 
to measure the relativistic advance of the periastron, which would likely be contaminated by 
classical effects as well. However, once a timing position is available, optical spectroscopy may 
allow determinarion of the mass ratio, as has been done in the case of PSR J 1740— 5340 in 
the globular cluster NGC 6397 ([271 . The feasibility of such a measurement will depend on the 
pulsar not being located in the dense and optically crowded core of the cluster. Fortunately, as 
PSR J 1748— 2446ad has a dispersion measure which is ~ 3 units lower than the ~ 239 pc cm -3 
average dispersion measure of the cluster pulsars, it is likely to be located outside of the core. 

Although there are selection effects, especially at radio wavelengths, against finding fast, bi- 
nary pulsars (21 ), we have maintained excellent sensitivity to these through the use of advanced 
search techniques (|22"|) and data with high time and frequency resolution. For example, these 
searches blindly detect the 596-Hz binary pulsar PSR J1748— 24460 (P or b ~ 6 hrs) at its second 
and fourth harmonic, which is equivalent to detecting a highly accelerating 1 192-Hz pulsar and 
its second harmonic. While these searches are clearly sensitive to pulsars faster than 716 Hz, 
we note that obscuration of the pulsar signal by material blown off the companion by the pulsar 
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wind may play an important role in reducing the chances of detecting such systems Of 
the five fastest known millisecond pulsars (including the two found in the Galactic plane; see 
Tab 2), four are eclipsing, and the fifth, PSR B 1937+21, is unusual because, unlike ~80% of 
Galactic plane millisecond pulsars, it has no binary companion. Since rotational energy loss is 
inversely proportional to the cube of the spin period (E = An 2 IP / P 3 , where E is the spin-down 
luminosity and P is the period derivative), it is plausible that a large fraction of the fastest bi- 
nary pulsars are evading detection because their powerful winds are ablating their companions. 
Some of the ablated material remains in the vicinity of the system and obscures radio pulsa- 
tions, particularly at lower radio frequencies. Although we have no detections of the pulsar at 
other frequencies, several other eclipsing pulsars are observed to have longer eclipse durations 
at lower radio frequencies If this is also the case for PSR J 1748— 2446ad, then the unusu- 
ally high observing frequency used here (2 GHz, while most globular cluster surveys have been 
conducted at 1.4 GHz or lower) was likely crucial in detecting this pulsar. PSR J1748— 2446ad 
is too weak (the flux density at 1950 MHz is ~0.08 mJy) to be detectable by the vast majority 
of Galactic plane surveys and its high eclipse fraction compounds this problem. This suggests 
that other even faster- spinning pulsars exist, but will require deeper surveys (perhaps at higher 
observing frequencies to mitigate obscuration of the pulsar signal by intra-binary material) and 
more concentrated efforts to be detected. In effect, the isolated nature and very large lumi- 
nosity of PSR B 1937+21 make it a unique object, rather than a representative member of the 
millisecond pulsar population. 

Low-mass X-ray binaries (LMXBs) with neutron star members are the likely progenitors 
of the radio millisecond pulsars. As the spin-up time for a neutron star to reach > 1000 Hz 
rotation via the accretion of matter in an LMXB is much shorter than the typical LMXB lifetime, 
one might naively expect many millisecond pulsars to be rotating at sub-millisecond periods. 
However, given the observed lack of pulsars spinning this fast, gravitational radiation has been 
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proposed as a limiting mechanism, as it could be responsible for carrying away rotational kinetic 
energy from the star, thus spinning it down. Specifically, gravitational waves may be acting 
either through an accretion-induced mass quadrupole on the crust (25), a large toroidal magnetic 
field (2<5j), or an r-mode (Rossby wave) instability in the stellar core (27 28). 

PSR J1748— 2446ad provides interesting constraints on the r-mode possibility. These os- 
cillations are believed to be present in all rotating neutron stars (29 30). Due to gravitational 
radiation emission, the r-modes become unstable and grow exponentially. The amplitude of 
the mode continues to grow as angular momentum is radiated away, and the star spins down. 
However, it is unclear whether the driving of these modes by gravitational waves can overcome 
viscous damping in the star. Damping depends on the core temperature of the star and its spin 
rate, as well as several other factors including the thickness of the neutron star crust and how 
it couples to the core. For a given core temperature, it is possible to derive a critical spin fre- 
quency above which the proposed mode is unstable and will cause the star to spin down rapidly. 
It has been predicted that the critical frequency is < 700 Hz for a wide range of core temperature 
(10 7 — 10 10 K) and a realistic model of the neutron star crust (28). Our discovery of a 7 16-Hz pul- 
sar indicates that if the r-mode instability limits neutron star spin-up, then either it must become 
important only at more rapid spin rates or better modelling of the neutron star crust is required. 
The current upper limit on the frequency derivative of PSR J 1748— 2446ad is consistent with 
those measured for other millisecond pulsars, and does not suggest an anomalously rapid spin- 
down rate. Thus, while there is currently no evidence that PSR J1748— 2446ad is spinning down 
due to gravitational radiation, the possible importance of systems like PSR J1748— 2446ad as 
gravitational wave sources for detectors like LIGO makes improved modelling of neutron star 
structure and gravitational wave instabilities critically important. 

The observed spin frequencies of LMXBs are all less than 620 Hz. The biases that exist 
at radio wavelengths against finding much faster pulsars do not exist for LMXBs, as X-rays 
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do not suffer from the dispersive effects of the interstellar medium and are less obscured by 
intra-binary material. This suggests that faster- spinning neutron stars in LMXBs should be 
detectable if they exist. However, transient coherent pulsations are only observed in seven 
sources, and there may be unidentified sources of bias in the population that are preventing faster 
pulsations from being detected. Chakrabarty et al. 2003 (4) performed a Bayesian statistical 
analysis on the spin frequencies of 1 1 nuclear-powered millisecond pulsars, those for which 
the spin frequency is known from the detection of burst oscillations, and concluded that the 
sample is consistent (at the 95% confidence level) with a cutoff v max = 760 Hz. More recent 
work by the same authors (U7T) . which added 2 pulsars to the sample, has revised this limit to 
730 Hz. Based on this, they conclude that something, possibly gravitational radiation, is limiting 
the spin frequency of neutron stars. If their statistically derived upper limit is realistic (and 
therefore some mechanism is limiting neutron star spin-up), then the 716-Hz pulsar presented 
here is likely an extremely rare object. However, we note that their Bayesian calculation is very 
sensitive to the neutron star with the slowest spin frequency included in the analysis. In addition, 
the assumption made by these authors that the pulsar spin rates are uniform in frequency (at 
least over some range vi ow to Uhigh) would likely not apply to the Terzan 5 pulsars, whose spin 
period distribution is clearly not uniform (1771) . even accounting for the observational bias against 
detecting the fastest pulsars. Hence, a recalculation of the maximum spin cutoff using the 
Terzan 5 sample of pulsars and the same statistical analysis would not be appropriate, although 
the existence of PSR J1748— 2446ad already implies that the cutoff must be higher. 
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Parameter 


Value 


Rotational Parameters 


Pulse period P (s) 


0.00139595482(6) 


Period derivative \P\ (s/s) 


< 6 x iir 19 


Pulse frequency v (Hz) 


716.35556(3) 


Frequency derivative z> (Hz/s) 


< 3 x 1(T 13 


Epoch (MJD) 


53500 


Orbital Parameters 


Orbital period P orb (days) 


1.09443034(6) 


Projected semi-major axis x (lt-s) 


1.10280(6) 


Time of ascending node Tasc (MJD) 


53318.995689(12) 


Eccentricity e 


< 0.0001 


Derived Quantities 


Companion minimum mass M 2)m in (M ) 


0.14 


Dispersion measure DM (pc cm~ 3 ) 


235.6(1) 


Flux density at 1950 MHz 51950 (mJy) 


0.08(2) 


Characteristic age r c (years) 


> 2.5 x 10 7 


Surface magnetic field B sur f (G) 


< 1.1 x 10 9 


Spin-down luminosity E erg/s 


< 1.3 x 10 37 



Table 1: Measured and derived parameters of PSR J1748— 2446ad. All measured spin and 
orbital parameters were determined using the TEMPO software package (U2"l). using arbitrary 
phase offsets between observing epochs. Given the currently sparsely sampled data, it is impos- 
sible to phase connect separate observations. For this reason, we provide only an upper limit on 
the magnitude of the spin frequency derivative of the pulsar, which incorporates the maximum 
possible acceleration due to the gravitational potential of Terzan 5 assuming a position close 
to the cluster center. Likewise, we can currently only place limits on the derived characteris- 
tic age, surface magnetic field and spin-down luminosity. The minimum companion mass is 
derived assuming a pulsar mass of 1.4 M Q . The dispersion measure was determined by measur- 
ing pulse arrival time delays in the different frequency channels across the 600-MHz observing 
bandwidth. 
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Spin Frequency 


P b 




Eclipse 




Pulsar 


(Hz) 


(days) 


(M ) 


Fraction 


Location 


J1748- 


-2446ad 


716.358 


1.0944 


0.14 


0.4 


Terzan 5 


B1937- 


m 


641.931 




isolated 




Galaxy 


B1957- 


f20 


622.123 


0.3819 


0.021 


0.1 


Galaxy 


J1748- 


-24460 


596.435 


0.2595 


0.035 


0.05 


Terzan 5 


J1748- 


-2446P 


578.496 


0.3626 


0.37 


0.4 


Terzan 5 


J1843- 


-1113 


541.812 




isolated 




Galaxy 


J0034- 


-0534 


532.714 


1.5892 


0.14 





Galaxy 


J1748- 


-2446Y 


488.243 


1.17 


0.14 





Terzan 5 


J1748- 


-2446V 


482.507 


0.5036 


0.12 





Terzan 5 


B0021- 


-72J 


476.048 


0.1206 


0.020 


0.1* 


47 Tucanae 



Table 2: The ten fastest- spinning known radio pulsars. Data compiled from the ATNF pulsar 
database (TJ?1). *B0021— 72 J is eclipsed only at radio frequencies <1 GHz. 
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0.5 1 

Pulse Phase 

Figure 1: Master PSR J1748— 2446ad pulse profile from the combination of eight GBT Pulsar 
Spigot observations at 2 GHz with particularly good detections of the pulsar. The cumulative 
integration time is ~54hrs. There are 32 phase bins across the profile, and the y-axis plots flux 
in arbitrary units. The one sigma error bar on the flux is shown in the lower left-hand corner. 
The effective time resolution of the data, ~ 300 /is, which accounts for pulse smearing due to 
channelisation of the dispersed data and finite time sampling, is indicated by the horizontal bar. 
A weak, but statistically significant interpulse is seen at a phase of ~0.75. 
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